Dynamics of momentum entanglement in lowest-order QED 
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We study the dynamics of momentum entanglement generated in the lowest-order QED interaction 
between two massive spin-1/2 charged particles, which grows in time as the two fermions exchange 
virtual photons. We observe that the degree of generated entanglement between interacting particles 
with initial well-defined momentum can be infinite. We explain this divergence in the context of 
entanglement theory for continuous variables, and show how to circumvent this apparent paradox. 
Finally, we discuss two different possibilities of transforming momentum into spin entanglement, 
through dynamical operations or through Lorentz boosts. 
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I. INTRODUCTION 

Entanglement is the physical property that 
Schrodinger described as Not one but the charac- 
teristic trait of quantum mechanics. It has played a 
fundamental role in the study of the completeness of 
quantum mechanics P, Q|. Nowadays, in quantum 
information theory, entanglement is considered as a 
physical resource, equivalent in many aspects to the 
role energy played in classical and quantum mechanics. 
Most applications in this novel field, like quantum 
teleportation quantum communication 4], quantum 
cryptography |5(, and some algorithms of quantum 
computation are carried out by using this intriguing 
quantum property @. A thorough study of entangle- 
ment in quantum information theory would demand a 
natural classification between discrete and continuous 
variables [8j. 

In the last few years two apparently different fields, 
entanglement and relativity, have experienced intense re- 
search in an effort for treat ing them in a common frame- 
work n in m 13 13 13 iii in m m. M 0S t f those 

works investigated the Lorentz covariance of entangle- 
ment through purely kinematic considerations, and only 
a few of them studied ab initio the entanglement dy- 
namics. For example, in the context of Quantum Elec- 
trodynamics (QED), Pachos and Solano [l4( considered 
the generation and degree of entanglement of spin cor- 
relations in the scattering process of a pair of massive 
spin-1/2 charged particles, for an initially pure product 
state, in the low-energy limit and to the lowest order in 
QED. Manoukian and Yongram 18] computed the effect 
of spin polarization on correlations in a similar model, 
but also for the case of two photons created after e + e~ 
annihilation, analyzing the violation of Bell's inequality 
0- In an earlier work, Grobe et al. studied, in the 
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nonrelativistic limit, the dynamics of entanglement in po- 
sition/momentum of two electrons which interact with 
each other and with a nucleus via a smoothed Coulomb 
potential. They found that the associated quantum cor- 
relations manifest a tendency to increase as a function of 
the interaction time. 

In this paper, we study to the lowest order in QED 
the interaction of a pair of identical, charged, massive 
spin-1/2 particles, and how this interaction increases the 
entanglement in the particle momenta as a function of 
time. We chose to work at lowest order, where en- 
tanglement already appears full-fledged, precisely for its 
simplicity. In particular this allows to set aside neatly 
other intricacies of QED, whose influence on entangle- 
ment should be subject of separate analysis. Here, the 
generation of entanglement is a consequence of a con- 
servation law: the total relativistic four-momentum is 
preserved in the system evolution. This will also be the 
case in any interaction verifying this conservation law, as 
occurs in closed multipartite systems, while allowing the 
change in the individual momentum of each component. 
In the asymptotic limit, the infinite spacetime intervals 
involved in the S-matrix result in the generation of an 
infinite amount of entanglement for interacting particles 
with well-defined momentum. QED is a place where in- 
finities can be avoided, and this will be also true, even 
though for other physical reasons, in the case of diver- 
gences appearing in momentum entanglement, a distinc- 
tive feature of continuous variables |g. We will also dis- 
cuss two different possibilities of establishing transfer of 
entanglement between momentum and spin degrees of 
freedom in the collective two-particle system: through 
dynamical operations or Lorentz boosts. 

In Sec. [n] we analize at lowest order and at finite time 
the generation of momentum entanglement between two 
electrons. In Sec. IIIII we calculate the Schmidt decom- 
position of the amplitude of a pair of spin-1/2 particles, 
showing the growth of momentum entanglement as they 
interact via QED. We obtain also analytic approxima- 
tions of the Schmidt modes, both, in momentum and 
configuration spaces. In Appendix ^ we include some 
notations and definitions related to entanglement theory 
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for discrete and continuous variables. In Appendix [B] 
we address the possibilities of transferring entanglement 
between momenta and spins via dynamical action, with 
Local Operations and Classical Communication (LOCC), 
or via kincmatical action, with Lorcntz transformations. 



II. TWO ELECTRON GREEN FUNCTION IN 
PERTURBATION THEORY 

To address the properties of entanglement of a two 
electron system one needs the amplitude (wave func- 
tion) ip(xi,x 2 ) of the system, an object with 16 spinor 
components dependent on the configuration space vari- 
ables Xi, X2 of both particles. The wave functions were 
studied perturbatively by Bethe and Salpeter |2fj | and 
their evolution equation was also given by Gell-Mann and 
Low |2l| . The wave function development is closely re- 
lated to the two particle Green function, 

K(l,2; 3, 4) = (^n^iM^MxaM^Mo) (1) 

which describes (in the Heisenberg picture) the sym- 
metrized probability amplitude for one electron to pro- 
ceed from the event X3 to the event x\ while the other 
proceeds from x^ to x 2 . If u ps (3) describes the electron 
at 3 and u p > s >(A) the one at 4, then 



iP(x 1 ,x 2 ) = I ^(3)^(4)^(1,25 3,4)7^7^ 



x u ps (3)itpv(4), 



(2) 



will be their correlated amplitude at 1, 2. In the free 
case this is just it ps (l) u p ' s i (2) , but the interaction will 
produce a reshuffling of momenta and spins that may 
lead to entanglement. The two body Green function K 
is precisely what we need for analysing the dynamical 
generation of entanglement between both electrons. 
Perturbatively [20|, 

Jr(l,2;3,4)=£jp(l,3)£fc(2,4) - e 2 j d A x 5 d i x 9 S F (l, 5) 

x S F (2, 6) tf B) D F (5, 6) 7(6)^(5, 3) S F (6, 4) 
+ ...-{1 <-> 2} (3) 

where all the objects appearing in the expansion are those 
of a free field theory. We may call to the successive 
terms on the right hand side of this expression. They will 
describe the transfer of properties between both particles 
due to the interaction. This reshuffling vanishes at lowest 
order, which gives just free propagation forward in time: 

d 3 x 1 d 3 x 2 d 3 x 3 d 3 x 4 u piS1 (1) u p2S2 (2)tf(°)(l, 2; 3, 4) 

X7°3)"p^(3)7(V Pb s 6 (4) = 9(ti-h)9(t 2 -U) 
x6 SlSa 5 S2Sb 6^( Pl - Po ) SW(p 2 - Pb ) (4) 

where u ps (x) = (2it)~ 3 t 2 (m/ E) 1 / 2 exp(— ipx)u s (p). The 
first effects of the interaction appear when putting 
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FIG. 1: Feynman diagrams for the QED interaction between 
two electrons (second order) . The minus sign denotes the an- 
tisymmetry of the amplitude associated to the fermion statis- 
tics. 



instead of K^ ' in the left hand side of the above equation. 
The corresponding process is shown in Fig. ^ To deal 
with this case we choose t\ = t 2 = t, t 3 = t± = — t and 
introduce the new variables 



(t 5 + t 6 ), t + e(-t,f). 



(5) 



*- = -t 6 ), t-e{-(t-\t+\),t-\t+\), (6) 

in Eq. ©, yielding 

(1,2; a,b;t) = 
2ie 2 



• u, ■ 
Jl a Jfj,2b 



(2tt) 4 y /2E 1 2E 2 2E a 2E b 
dk° 



0)2 
t-\t+\ 



^(fc ) 2 - (p„ - Pi) 



(5 (3) (P1+P2-Pa-Pb) 



dt -HE a +E b - El -E 2 )t + 



ie 



-t 



< / dt _ e -i{E,-E 2+ E b -E a +2k°)t- _ {1 ^ 2} , (7) 

>-(t-\t+\) 

where K^(l, 2; a, b; t) is a shorthand notation for what 
corresponds to (0J at second order, and j'jy = UkJ^ui, 
After some straightforward calculations, we obtain 

£00(1, 2; a, 6; t) = * ^ ^ + *> 2 Z P f^ l 
v y 4tt 3 y/2E 1 2E 2 2E a 2E b 

xULiMSt® + r t (t)]-jSj„ lb [s v (t) + r u (t)}}, 

(8) 



with, 
St(t) = 

Tt(*) = 



: EM 2 Eh - Ea ) -(p a -pi) 2 

sin [(Eh. +E 2 -E a - E b )t] 



Ei + E 2 — E a — E b 



1 



Pa - Pi I 



X < I 



1 



1 



^(E 2 -^ 2 ) v(Y?-v 2 ) 
1 1 

cos(Ei) 



E 2 — [i 2 Yl 2 -v 2 
1 



E 2 — [i 2 



-ifit 



1 



E 2 



(9) 



E sin(Ei) 



, (10) 
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/' 
v 
A 



and 



S u (t) 



Ei + E 2 — E a — Eb, 
A + 2| Po -pi|, 
-A + 2|p Q - Pl |, 

Ei — Ei + Eb — E a , 



S t (t),T tt (t)«->T t (t), 
1 2 



(11) 
(12) 

(13) 
(14) 



(15) 



St,u are the only contributions that remain asymptoti- 
cally (t — > oo) leading to the standard scattering am- 
plitude, while T t;U vanish in this limit. We recall that 
these are weak limits: no matter how large its modulus, 
the expression in Eq. (|10fl will vanish weakly due to its 
fast oscillatory behavior. On the other hand, the sine 
function in Eq. JUJl enforces energy conservation via 



lim 



sin[(Ei +E 2 -E a -E b )t] 



*oo Ei + E2 — E a — Eb 



= ir5(E 1 +E 2 -E a -E b 



(ie) 

This limit shows also that the entanglement in energies 
increases with time , see Appendix reaching its 
maximum (infinite) value when t — > 00 for particles with 
initial well-defined momenta and energy. This result is in- 
dependent of the chosen scattering configuration. Exact 
energy conservation at large times, united to a sharp mo- 
mentum distribution of the initial states, would naturally 
result into a high degree of entanglement. The better de- 
fined the initial momentum of each electron, the larger 
the asymptotic entanglement. The physical explanation 
to this unbounded growth is the following: the parti- 
cles with well defined momentum (unphysical states) are 
spread over all space, and thus their interaction is ubiqui- 
tous, with the consequent unbounded degree of generated 
entanglement. This is valid for every experimental setup, 
except for those pathological cases where the amplitude 
cancels out, due to some symmetry. In the following sec- 
tion, and for illustrative purposes, we will single out these 
two possibilities. 

i) The case of an unbounded degree of attainable entan- 
glement due to an incident electron with well defined mo- 
mentum. We consider, with no loss of generality, a fuzzy 
distribution in momentum of the second initial electron, 
for simplicity purposes. 

ii) Basically the same setup as in (i) but with a spe- 
cific spin configuration, which leads to cancellation of the 
amplitude at large times due to symmetry, and thus to 
no asymptotic entanglement generation. 

On the other hand, for finite times, nothing prevents 
a sizeable contribution from Eq. (|10fl . In fact, in the 
limiting case where i -1 is large compared to the ener- 
gies relevant in the problem, it may give the dominant 
contribution to entanglement. Whether the contribution 
from Tt(t) and T u (t) is relevant, or not, depends on the 
particular case considered. 



III. TWO ELECTRON ENTANGLEMENT 
GENERATION AT LOWEST ORDER 

The electrons at 13 , X4 will be generically described by 
an amplitude F 

ip F (x 3 ,X4) = ^2 rf3pa d3 PbF{p ai s a ;p bl s b ) 

xu p a ,Sa( x 3)u pb , Sb (x4) (17) 

that should be normalizable to allow for a physical inter- 
pretation, i.e., 



E 



d 3 p a I d J p fc |F(p a , s a ; p bl s b ) 



1. 



(18) 



For separable states where F(a; b) = f a (p a , s a )f b (p b , s b ), 
f a and fb could be Gaussian amplitudes g centered 
around a certain fixed momentum p° and a certain spin 
component s , 



-(p- 



5 a) 3 / 2 



which in the limit of vanishing widths give the standard 
-well defined- momentum state S ss o8^ 3 \p — p°). 

In the absence of interactions, a separable initial state 
will continue to be separable forever. However, interac- 
tions destroy this simple picture due to the effect of U. 
Clearly, the final state 

F {2) (pi, si; p 2 , s 2 ; t)=Y^ J d3 P*J d 3 PbK (2) (1, 2; a, b; t) 



x F(p a ,s a ;Pb,Sb) 



(19) 



can not be factorized. 

In the rest of this section we analyze the final state 
F^ 2 ' (pi, si; P2, S2', t) in Eq. 1)19(1 to show how the vari- 
ables pi and P2 get entangled by the interaction. We 
consider the nonrelativistic regime in which all interven- 
ing momenta and widths p, a -C m, so the characteristic 
times t under consideration are appreciable. We single 
out the particular case of a projectile fermion a scattered 
off a fuzzy target fermion 6 centered around p£ = 0. As a 
further simplification, we consider the projectile momen- 
tum sharply distributed around p" (a a <C p„) so that 
the initial state can be approximated by 



F(a;b)^6 SaS o6^(p a -p° a )- 



(v/f^) 3/2 



-(Pb 



Our kinematical configuration would acquire complete 
generality should we introduce a finite momenta p° for 
the initial electron b. The reference system would be in 
this case midway between the lab. system and the c.o.m. 
system. In short, the choice p^ = will not affect the 
qualitative properties of entanglement generation. 
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FIG. 2: Experimental setup considered in the calculations. 



where £ 



£(p + q), p, = v^P 2 + (* 



and 



< = Y^-t. In the following, we analyze different spe- 
cific spin configurations in the non-relativistic limit with 
the help of Eq. i|23|l . We consider an incident particle 
energy of around 1 eV<§; m (p Q — 1 KeV), and a mo- 
mentum spreading a one order of magnitude less than 
p\. We make this choice of p° a and a to obtain longer 
interaction times, of femtoseconds (t = pj^t). Thus the 
parameter values we consider in the subsequent analysis 
are p°/m = 0.002 and a/m = 0.0002. We consider the 
initial spin state for particles a and b as 



(24) 



We will work in the lab frame, where particle 6 shows a 
fuzzy momentum distribution around = 0, and focus 
in the kinematical situation in which the final state mo- 
menta satisfy Pi • P2 = and also p a ■ p" = l/\/2p a p®, 
a = 1 , 2 (see Fig. |2} . This choice not only avoids forward 
scattering divergencies but also simplifies the expression 
of the amplitude in Eq. (|19fl . due to the chosen angles. 
For sure, the qualitative conclusions would also hold in 
other frames, like the center-of-mass one. We obtain 



Ei + E 2 — E a 



E, 



|P6=Pl+P2 — P a 



P°a 



(Pl + P2 



p 2 -V2p a ) + O((p a /m) 3 p°a), 
-p°) 2 ( Pl -pl/V2) 2 , (pa-/,. 



(Pl - Pa) 2 = (Pl "PS 
(P3 - Pa)' = (P2 - Pa 



+ (Pa) 2 /2, 
2 + (Pa) 2 /2. 



(21) 



Here, boldface characters represent trivectors, otherwise 
they represent their associated norms. We perform now 
the following change of variables, 



p 1 

T2 = * {! " 



V2 



g 

V2 



- 1 

(7 



,22) 



along an arbitrary direction that will serve to mea- 
sure spin components in all the calculation. The phys- 
ical results we are interested in do not depend on 
this choice of direction. The QED interaction, in the 
non-relativistic regime considered, at lowest order, is 
a Coulomb interaction that does not change the spins 
of the fermions. In fact, {jfj^b) ^ 4m 2 i s » Sl i5 s o S2 , 
(jibjfi2a) — 4m # s o sl <5 s o S2 . Given the initial spin states 
of Eq. I|24fl . depending on whether the channel is t or u, 
the possible final spin states are 



\SlS 2 )t 

\sis 2 )u 



IT). 



(25) 
(26) 



Due to the fact that the considered fermions are identical, 
the resulting amplitude after applying the Schmidt proce- 
dure is a superposition of Slater determinants [23 . 124L |25I | . 
Whenever this decomposition contains just one Slater de- 
terminant (Slater number equal to 1) the state is not en- 
tangled: its correlations are just due to the statistics and 
are not useful for the applications because they do not 
contain any additional physical information. If the am- 
plitude contains more than one determinant, the state is 
entangled. Splitting the amplitude in the corresponding 
ones for the t and u channels, we have 



turning the amplitude in Eq. (|19(l into 



F {2) {p,si;q,s 2 ;t) oc 



sin[(p + q)t] 



4- 



(jlaJ^b) Sb=s S (jii)jfj,2a) s 



A/2 



1 



V/2 e -o 2 /2 



£ 2 



cos[(p + q)t] - e 



A 2 ) 



tsm[(p + q)i\ 



- fit 



{p,l ^q,2} e- p2 / 2 e"« 2 / 2 , 



(23) 



with 



F( 2 )(p,T; g ,|;t) t cx 
sin[(p 4- q)i\ 1 



- P V2 e -oV2 



M/2 I /i(£2 _ ^2) 



p2 + (^y 
1 

T,sm[(p + q)t] 



£ 2 -ft 2 



cos[(p 4- - e ^ 



x e -P 2 /2 e -« 2 /2 f 



(27) 



FW( P ,l;q,V,t) u ^ FW(p,v,q,l;t) t , 

p^q. (28) 
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In the infinite time limit the sine function converges to 
5(p + q), which is a distribution with infinite entangle- 
ment [23 ■ The presence of the sine function is due to the 
finite time interval of integration in Eq. Q. This kind 
of behavior can be interpreted as a time diffraction phe- 
nomenon [2l| . It has direct analogy with the diffraction 
of electromagnetic waves that go through a single slit of 
width 2L comparable to the wavelength A. The analogy 
is complete if one identifies t with L and p + q with 2tt/X. 

In Fig. |3 we plot the modulus of Eq. (|2"7|) versus p, q, 
at times t = 1, 2, 3, 4. This graphic shows the progressive 
clustering of the amplitude around the curve q = —p, 
due to the function S1SK£±£M_ This is a clear signal 

p+q b 

of the growth in time of the momentum entanglement. 
Fig. |3 puts also in evidence the previously mentioned 
time diffraction pattern. 



t 



t = 1 



o 



3 

t = 4 




2 p 2 
t = 3 



t, p z 



FIG. 3: (Color online) |F (2) (p, |; q, I; t)t\ versus p, q at t = 
1,2,3,4 . 

We have applied the method for obtaining the Schmidt 
decomposition given in Ref. j22J to Eq. (|27[) . considering 
for the orthonormal functions {0^\p)}, {0^ 2 \q)} Hcr- 
mite polynomials with their weights, to take advantage 
of the two Gaussian functions. We obtain the Schmidt 
decomposition for t = 1,2,3,4, where the error with ma- 
trices C mn 12 x 12 or smaller is d^ n ng < 7 ■ 10~ 3 in all 
considered cases. We plot in Fig. [tj the coefficients A„ 
of the Schmidt decomposition of Eq. (|27|l as a function 
of n, for times t = 1,2,3,4. The number of A„ differ- 
ent from zero increases as time is elapsed, and thus the 
entanglement grows. 

The complete Schmidt decomposition, including chan- 
nels t and u, is given in terms of Slater determinants |23| . 



12345012345 

n n 



t = 3 



t = 4 




FIG. 4: Eigenvalues X n versus n at times t = 1,2, 3, 4. 



and is usually called Slater decomposition. It is obtained 
antisymmetrizing the amplitude for channel t 



F^(p, Sl ;q,s 2 ;t)^Y,\/ X ^ 



^\p,i)\ TW(9,*)| l)-WM\ l)^>{q,i)\ t) 



i.(2) 



(2), 



.,(1), 



V2 



(29) 



where the modes ip„ (k, t) and ipn^ (k, t) are the Schmidt 
modes of the channel t obtained for particles 1 and 2 
respectively, and they correspond to the modes of the 
channel u for particles 2 and 1 respectively. 

A measure of the entanglement of a pure bipartite state 
of the form of Eq. I|29|l. equivalent to the entropy of en- 
tanglement S, is given by the Slater number 



ropj 

m 



K 



1 



(30) 



K gives the number of effective Slater determinants 
which appear in a certain pure bipartite state in the form 
of Eq. (|29l) . The larger the value of K, the larger the en- 
tanglement. For K = 1 (one Slater determinant) there 
is no entanglement. This measure is obtained consider- 
ing the average probability, which is given by X^^Lo X n 
E^Lo ^™ = 1> an d tnus {^n} can be seen as a probabil- 
ity distribution) . The inverse of the average probability 
is the Slater number. Its attractive properties are that it 
is independent of the representation of the wavefunction, 
it is gauge invariant, and it reaches its minimum value 
of 1 for the separable state (single Slater determinant). 
In Fig. we show the Slater number K as a function of 
elapsed time t, verifying that the entanglement increases 
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n = 



n = 1 




-0.8 




n = 2 



n = 3 



FIG. 5: Slater number K as a function of the elapsed time t. 



as the system evolves. It can be appreciated in this fig- 
ure the monotonic growth of entanglement, due to the 
fact that we have considered an incident electron with 
well defined momentum. In realistic physical situations 
with wave packets, this growth would stop, due to the 
momentum spread of the initial electrons. The general 
trend is that the higher the precision in the incident elec- 
tron momentum, the larger the resulting asymptotic en- 
tanglement. The fact that the entanglement in momenta 
between the two fermions increases with time is a conse- 
quence of the interaction between them. We remark that 
the entanglement cannot grow unless the two particles 
"feel" each other. The correlations in momenta are not 
specific of QED: the effect of any interaction producing 
momentum exchange while conserving total momentum 
will translate into momentum correlations. 

The Schmidt modes in momenta space for the ampli- 
tude of Eq. l(77|) are given by 

"0 



n=0 



a = 1,2, 



(31) 



where uq is the corresponding cut-off and the values of 
the coefficients im» (t) are obtained through the method 
given in Ref. |22j . The modes in momenta space de- 
pend on time because they are not stationary states: the 
QED dynamics between the two fermions and the inde- 
terminacy on the energy at early stages of the interaction 
give this dependence. By construction, the coefficients 
Amn(t) do not depend on p, q. 

We plot in Fig.|51the Schmidt modes ijfflip, t) at times 
t= 1, 2, 3, 4 for n = 0, 1, 2, 3 (we are plotting specifically 
the real part of each mode only, which approximates well 
the whole mode, because Eq. (I27|l is almost real for the 
cases considered). The sharper modes for each n cor- 
respond to the later times. Each Schmidt mode is well 
approximated at early times by the corresponding Her- 
mite orthonormal function, and afterwards it sharpens 
and deviates from that function: it gets corrections from 
higher order polynomials. The fact that the modes get 
thinner with time is related to the behavior of Eq. I|27|) 
at large times. In particular the sine function goes to 




FIG. 6: (Color online) Schmidt modes ipn\p,t) at times 
t = 1,2,3,4 for n = 0,1,2,3. The sharper modes for each 
n correspond to the later times. 



5(p + q) and thus the amplitude gets sharper. 

Now we consider the Schmidt modes in configuration 
space. To obtain them, we just Fourier transform the 
modes of Eq. Q31f) with respect to the momenta p\ , p2 



^\x a ,t) = 



1 



d Po e*fe-«-a**tyW(A(Pa),t), 

where a — 1,2. The dependence of p on pi and q on p 2 is 

given through Eq. lj2*2|) . The factor e -1 ^' in Eq. (J32J is 
the one which commutes the states between the interac- 
tion picture (considered in Eq. (|31|l and in the previous 
calculations in Sees. [H]and[ffi]) and the Schrodinger pic- 
ture. 

The Hermite polynomials obey the following expres- 
sion El 



dxe- {x - y)2 H n {ax) = V^(l-a 2 ) n/2 H n 



ay 



vl — a 

(33) 

With the help of Eq. i|33|) and by linearity of the Fourier 
transforms, we are able to obtain analytic expressions 
for the Schmidt modes in configuration space (to a cer- 
tain accuracy, which depends on the cut-offs considered) . 
This is possible because the dispersion relation of the 
massive fermions in the considered non-relativistic limit 
is E a — and thus the integral of Eq. (|32() can be 
obtained analytically using Eq. I|33|) . 

The corresponding Schmidt modes in configuration 
space are then given by 

^\x a ,i) ~ A^ n (t)Oi a \x a ,t), a = 1,2, (34) 
n=0 
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where the orthonormal functions in configuration space 
are 



-in arctan(crf )+-icr (x a —t/2) 



(35) 



In Eqs. (|3~H) and we are using dimensionless vari- 

ables, x a — a = 1,2, a = p-, and the dimensionless 

time defined before, t = %m*' The modes in Eqs. i|34|) 
and l|35fl are normalized in the variables x a . The or- 
thonormal functions of Eq. Q35[) propagate in space at 

a speed Jr£ m and they spread in their evolution. Addi- 
tionally, the modes of Eq. (|34|) have also the time depen- 
dence of Ami, (t) . The Slater decomposition in configura- 
tion space, obtained Fourier transforming the modes of 
Eq. EH> is 



[^(£i,t)| T)^ a; (5a,t)l 4) 



^i^)l l>^ 2 J)l T}]- 



(36) 



The coefficients X n (t) are unaffected by the Fourier trans- 
formation, and thus the degree of entanglement in con- 
figuration space is the same as in momenta space. 
We consider now the initial spin configuration 



\s°A) = | TT), 



(37) 



where, the only possible final state in the non-relativistic 
limit is 



\sis 2 ) = | TT)- 



(38) 



In this case, the sine term goes to zero, because the mo- 
mentum part of this term is antisymmetric in p 2 , q 2 and 
the sine function goes to 5{p + q), which has support (as 
a distribution) on q = —p. We point out that the sine 
contribution to this amplitude is negligible because of the 
particular setup chosen. In other experiment configura- 
tions the amplitude in Eq. (|19fl associated to the spin 
states of Eqs. (|37|) and (|38|l would have appreciable sine 
term and thus increasing momenta entanglement with 
time. On the other hand, in this case the contribution 
from T t (t) in Eq. (|1L)|) and T u (t) is even smaller than the 
sine term, and converges weakly to zero. 

We plot in Fig. \7\ the real and imaginary parts of the 
term associated to T t (t) and T u (t) in Eq. I|23[l . which we 
denote by g(p, q, t), for spin states of Eqs. i|57|) and 
as a function of time t <G (1,1.001) and having p = 1, 
q = 1.2. We want to show with it the strong oscilla- 
tory character of the amplitude with time, and how all 
the contributions interfere destructively with each other 
giving a zero final value. This is similar to the station- 
ary phase procedure, in which only the contributions in 



proximity to the stationary value of the phase do interfere 
constructively and are appreciable. What we display here 
is the weak convergence to zero for the functions T t (i) 
and T u (t). 




= 1.2, t)] 



FIG. 7: Real and imaginary parts of the amplitude g(p,q,t) 
for p = 1, q = 1.2, and t 6 (1, 1.001) (arbitrary units). 

In this section, we investigated the generation of en- 
tanglement in momenta between two identical spin-1/2 
particles which interact via QED. We showed how the 
correlations grow as the energy conservation is increas- 
ingly fulfilled with time. The previous calculation had, 
however, the approximation of considering a projectile 
particle with perfectly well defined momentum, some- 
thing not achievable in practice. This is a first step to- 
wards a real experiment, where both fermions will have 
a dispersion in momenta and thus infinite entanglement 
will never be reached, due to the additional integrals of 
the Dirac delta S(AE) over the momentum spread. 



IV. CONCLUSIONS 

We analyzed the dynamical generation of entanglement 
between two electrons due to their mutual interaction by 
means of the Slater number. In the asymptotic limit, and 
for electrons whose initial momenta are sharply defined, 
entanglement divergencies may appear. We observe that 
considering finite-time intervals of interaction, and/or a 
certain spreading in the particles momentum, the entan- 
glement remains finite. We have studied for the first time 
the dynamical generation of momentum entanglement of 
two spin-1/2 particles at lowest order QED, observing 
how the correlations increase as the particles exchange 
virtual photons. We obtain the Schmidt decomposition 
of the scattering amplitude, given in terms of Slater de- 
terminants, and the Slater number, which clearly shows 
the growth of entanglement with time. We also obtain 
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analytic approximations of the Schmidt modes, both in 
momentum and configuration spaces. 
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APPENDIX A: ENTANGLEMENT IN FINITE 
AND INFINITE DIMENSIONAL HILBERT 
SPACES 

We consider a composite system S described by a 
Hilbert space H, which may be either finite or infinite 
dimensional. This space is constructed as the tensor 
product of the Hilbert spaces associated to each of the 
subsystems, S a , of S. For simplicity, and in view of our 
present purposes, we restrict ourselves to pure states of 
a bipartite system S. Thus, a = 1, 2, and H = Hi <£> H 2 - 

Definition: product state. A vector state |\&) of the 
system S is a product state if it can be written as 

|*) = |*( 1 ))|*W), (Al) 

where |*C0) e H x and |^ 2 )) E H 2 - 

Definition: entangled state. A vector state \^S) of the 
system S is entangled if it is not a product state. A 
relevant example, where dim(7ii) = dim(7i2) = 2, is 
called the singlet state, 

\^-) = ^\^ ) )\^ ) )-H 1) )\^ ) )). (A2) 

A very useful tool for analyzing the entanglement of pure 
states of bipartite systems is given by the Schmidt decom- 
position [28L |2^| . It basically consists in expressing the 
pure bipartite state as a sum of biorthonormal products, 
with positive coefficients VAn, as follows 

d-i 

l*)=Ev / ^l*l 1) >l*i 2) ) I (A3) 
n=0 

where {|\E'! 1 ' ) )}, {l^n^)}, are orthonormal bases associ- 
ated to Hi and H2, respectively. In Eq. (|A3(I . d — 
min{dim(7ii), dim(7i2)}, and it may be infinite, as for 
in continuous variable systems, describing momentum, 
energy, position, frequency, or the like. In those cases, 
the states ) would be (square integrable) L 2 wave 
functions, 

(p|*W) = ^(p), a = 1,2, (A4) 



where p denotes the corresponding continuous variable. 

For pure bipartite states a relevant measure of entan- 
glement is the entropy of entanglement, S. Given a state 
\^>), it is defined as the von Neumann entropy of the 
reduced density matrix with respect to Si or S 2 , 

S = -^2 A„log 2 A„, (A5) 

n=0 

where the A„'s are given in Eq. (|A3() . In general, S > 0, 
5 = for a product state, and the more entangled a 
state is, the larger S. For a maximally entangled state, 
S = log 2 d, and if d = 00, then S diverges. 

An interesting work where the Schmidt decomposition 
for continuous variables is analyzed discretizing the cor- 
responding integral equations can be found in Ref. |29j |. 
Another method for obtaining the continuous variables 
Schmidt decomposition, based in decomposing the bi- 
partite wave function in complete sets of orthonormal 
functions, is developed in Ref. p2^ . 



APPENDIX B: ENTANGLEMENT TRANSFER 
BETWEEN MOMENTUM AND SPIN 

1. Dynamical transfer and distillation 

In Sec lIIII we computed the entanglement in momenta 
for a pair of identical spin- 1/2 particles which interact 
through exchange of a virtual photon. The sharper the 
initial momentum distribution of the incident fermion, 
and the longer the interaction time, the larger the entan- 
glement in momenta. Heisenberg's principle, on the other 
hand, establishes a limit to the precision with which the 
momentum may be defined and hence to the achievable 
degree of entanglement. 

It is possible, in principle, to transform the entangle- 
ment in momenta into entanglement in spins. This is 
easily seen in terms of the majorization criterion |rll30|. 
which is of practical interest because the experimental- 
ist usually manipulates spins. Here, we will analyze this 
entanglement transfer. 

Majorization is an area of mathematics which pre- 
dates quantum mechanics. Quoting Nielsen and Chuang, 
"Majorization is an ordering on d-dimensional real vec- 
tors intended to capture the notion that one vector is 
more or less disordered than another". We consider 
a pair of d-dimensional vectors, x = (x\, X4) and 
y = (yi, yd). We say x is majorized by y, written 
x-<y,tf Ej=i x j ^ Ej=i Vj for k = l, d, with equal- 
ity instead of inequality for k = d. We denote by z* the 
components of z in decreasing order (zj > z\ > ... > z^). 
The interest of this work in the majorization concept 
comes from a theorem which states that a bipartite pure 
state may be transformed to another pure state \<f>) by 
Local Operations and Classical Communication (LOCC) 
if and only if A^ -< A^, where A^,, A^ are the vectors 
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of (square) coefficients of the Schmidt decomposition of 
the states \<fr), respectively. LOCC adds to those 
quantum operations effected only locally the possibility 
of classical communication between spatially separated 
parts of the system. According to this criterion, it would 
be possible in principle to obtain a singlet spin state \(f>) 
beginning with a momentum entangled state \ip) when- 
ever Xjf, -< A^. 

The possibility of obtaining a singlet spin state from 
a momentum-entangled state can be extended to a more 
efficient situation: the possibility of distillation of entan- 
glement. This idea consists on obtaining multiple singlet 
states beginning with several copies of a given pure state 
\t/j). The distillable entanglement of \ip) consists in the 
ratio n/m, where m is the number of copies of we 
have initially, and n the number of singlet states we are 
able to obtain via LOCC acting on these copies. It can 
be shown p that for pure states the distillable entan- 
glement equals the entropy of entanglement, S. Thus, in 
the continuous case (infinite-dimensional Hilbert space), 
the distillable entanglement is not bounded from above, 
because neither is S. According to this, the larger the en- 
tanglement in momenta the more singlet states could be 
obtained with LOCC. 

To illustrate the possibility of entanglement trans- 
fer with a specific example, we consider a momentum- 
entangled state for two distinguishable fermions 

\i>) = -j=[4 x) (p)4 2 Hi) +4 1) (p)4 2 \i)}\ tt). (bi) 

This state has associated a vector A^, = 
(1/2,1/2,0,0,...). On the other hand, the singlet 
state 

l0>=4 1) Wi 2) (g)^(ITI>-|IT» (B2) 

has associated a vector A^ = (1/2, 1/2, 0, 0...). 

These vectors obey A^, -< A^, and thus the state en- 
tangled in momenta may be transformed into the state 
entangled in spins via LOCC. 

2. Kinematical transfer and Lorentz boosts 

Another approach to the study of entanglement trans- 
fer between momentum and spin degrees of freedom is 
the kinematical one. In fact, the Lorentz transforma- 
tions may entangle the spin and momentum degrees of 
freedom. To be more explicit, and following the notation 
of Ref. [l2T |. we consider a certain bipartite pure wave 
function g\ a (p, q) for two spin-1/2 fermions, where A and 
a denote respectively the spin degrees of freedom of each 
of the two fermions, and p and q the corresponding mo- 
menta. This would appear to an observer in a Lorentz 



transformed frame as 

5A.( P ,q) ^ ^r^^WtA-^A-iq), 
X'u' 

(B3) 

where 

cr$ = D$ 2 >(ie(A,p)) (B4) 

is the spin 1/2 representation of the Wigner rotation 
i?(A,p). The Wigner rotations of Eq. (|B4J) can be seen 
as conditional logical operators, which rotate the spin a 
certain angle depending on the value of the momentum. 
Thus, a Lorentz transformation will modify in general 
the entanglement between momentum and spin of each 
individual electron. We distinguish the following three 
cases. 

i) Product state in all variables. In this case, 

9Xa(p,q) =#i(p)p 2 (q)|A)|o-), (B5) 

and the entanglement at the rest reference frame is zero. 
Under a boost, the Wigner rotations of Eq. (|B4|) entangle 
the momentum of each fermion with its spin, and thus 
the entanglement momentum-spin grows |10| . 

ii) Entangled state spin-spin and/or momentum- 
momentum. We consider now a state 

SA CT (p,q) = /(p,q)|0) (B6) 

with |0) an arbitrary state of the spins, and /(p, q) 
an arbitrary state of the momenta. In this case, a 
Lorentz boost will entangle in general each spin with its 
corresponding momentum, and a careful analysis shows 
that the spin-spin entanglement never grows |l2j |. Of 
course, by applying the reversed boost the entanglement 
momentum-spin would be transferred back to the spin- 
spin one, and thus the latter would grow. This particular 
case shows that, for the state we considered in Sec. 11111 
given by Eqs. (|24|) . I|25|l and l|27|l . the entanglement could 
not be transferred from momenta into spins via Lorentz 
transformations. Thus, the dynamical approach would 
be here more suitable. 

in) Entangled state momentum- spin. According to 
the previous theorem, the momentum-spin entanglement 
may be lowered, transferring part of the correlations to 
the spins, or increased, taking some part of the correla- 
tions from them. To our knowledge, there is not a similar 
result for the momentum, that is, whether the momen- 
tum entanglement can be preserved under boosts, or it 
suffers decoherence similarly to the spins, and part of it 
is transferred to the momentum-spin part. This is a very 
interesting question, which we will treat more deeply in 
future works. 
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